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1  Introduction 


With  the  recent  advances  in  fast  algorithms  (see,  for  example,  [2])  for  solving  integral  equations 
resulting  from  potential  theory,  there  is  a  renewed  interest  in  the  classical  potential  theory.  As 
is  well  known,  the  jump  relations  of  the  single  and  double  layer  potentials  play  an  important 
role  in  the  classical  potential  theory.  In  [3],  the  jump  relations  of  the  quadruple  layer  potential 
on  a  curve  in  two  dimensions  are  derived;  and  it  is  shown  that  the  jumps  are  proportional  to  the 
product  of  the  density  of  the  potential  and  the  curvature  of  the  curve.  In  this  note,  we  derive 
the  jump  relations  of  the  quadruple  layer  potential  on  a  regular  surface  in  three  dimensions 
and  show  that  the  jumps  of  the  quadruple  layer  potential  are  proportional  to  the  product  of 
the  density  and  the  mean  curvature  of  the  underlying  surface.  The  result  is  summarized  in 
Theorem  3.8. 

2  Analytical  Preliminaries 

In  this  section,  we  collect  some  well  known  facts  from  classic  analysis  to  be  used  in  the  remainder 
of  the  paper. 

2.1  Notation 

We  will  denote  by  5  a  sufficiently  smooth  (say,  at  least  twice  continuously  differentiable)  regular 
surface  in  R3.  When  S  is  an  open  surface  with  boundary,  we  assume  that  S  is  an  open  set,  i.e., 
S  does  not  contain  its  boundary  C.  For  a  point  a;  €  S',  we  denote  the  unit  normal  vector  to  S 
at  x  by  N(x).  For  a  vector  x  =  (xi,X2,xs)  €  R3  we  will  denote  its  length  by  |x|.  Finally,  for 
two  vectors  x,y  €  M3,  we  denote  their  inner  product  by  (x,  y). 

2.2  Single,  Double,  and  Quadruple  Layer  Potentials 

For  any  x,  t  £  R3  and  x  ^  t,  we  define  the  Green’s  function  for  the  Laplace  equation  in  R3  via 
the  formula 


1 


Suppose  now  that  S  is  a  sufficiently  smooth  regular  surface.  For  t  e  S  we  consider  the  direc¬ 
tional  derivatives  of  the  function  G  with  respect  to  t  along  the  normal  directions  of  S  at  t.  It 
is  easy  to  verify  that 


9G(*'()-Ar(t).V«C0M)= 


dN(t ) 
d2G(x,t) 


x  ■ 


tl3 


=  N(t)-VtVtG(x,t)-N(t) 


3  {N(t),x  —  t)‘ 


(2) 


(3) 


9N(t) 2  |ar_t|5  |a;_t|3- 

In  the  hterature,  (1),  (2),  (3)  are  often  referred  to  as  the  single,  double,  and  quadruple  potentials 
respectively.  Suppose  further  that  a  :  S  M  is  a  sufficiently  smooth  function.  We  will  refer 
to  the  functions  given  by  the  formulae 


/  G(x,  t)  •  aft)  ■  dt , 

(4) 

Js 

f  dG(x,  t)  ]x 

Js  dN(t)  '  '  di' 

(5) 

f  d2G(x,t)  „ 

Js  dN(t) 2  (  )  dt ’ 

(6) 

as  the  single,  double,  and  quadruple  layer  potentials,  respectively. 


2.3  Finite  Part  Integrals  on  a  Regular  Surface  in  R3 

Finite  part  (also  referred  to  as  hypersingular)  integrals  on  Euclidean  spaces  are  well  known  and 
extensively  used  in  mechanical  engineering.  Here  we  generalize  the  definition  for  the  Euclidean 
spaces  given  by  Samko  (see  Chapter  1,  Section  5.3  in  [4])  to  define  finite  part  integrals  on  a 
regular  surface  in  R3 . 


Definition  2.1  Suppose  that  S  is  a  sufficiently  smooth  regular  surface.  Suppose  further  that 
the  function  f  is  integrable  in  S  -  De(x)  for  all  e  €  (0,  e0)  for  some  eo  >  0,  where  Dt( x)  =  {t  G 
S'Ht  —  x\  <  e}.  Then  f  is  said  to  possess  the  Hadamard  property  at  x  if  there  exist  constants 
afc,  b,  and  A*,,  possibly  complex-valued,  but  Re(Xk)  >0,  k  =  1, 2, . . . ,  N,  such  that 


L 


S—Dc(x) 


N  1 

f{t)dt  =  ak  ■  e  Xk  +  b-  log  -  +  J0(e), 

fc=  l  e 


(7) 
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where  lime_+o  /o(e)  exists  and  is  finite.  In  this  case  the  finite  part  of  the  (divergent)  integral 
Is—De  (x)  f(t)dt  is  defined  by  the  formula 

f-P-  [  f(t)dt  =  limlo(e).  (8) 

Js  e~*° 

Remark  2.1  When  S'  is  a  flat  surface  embedded  in  R2,  the  above  definition  coincides  with  the 
conventional  definition  for  finite  part  integrals  in  Rn  given  in  [4]. 

2.4  Local  Properties  of  a  Regular  Surface  in  R3 

It  is  well  known  that  locally  any  sufficiently  smooth  regular  surface  in  R3  is  the  graph  of 
a  sufficiently  smooth  function.  When  the  local  coordinates  are  chosen  along  the  principal 
directions,  the  surface  admits  a  particularly  simple  parametrization.  The  following  lemma 
summarizes  this  fact;  it  can  be  found  in  [1], 

Lemma  2.1  Suppose  that  S  is  a  sufficiently  smooth  regular  surface  in  R3  and  that  x  is  a  point 
in  S.  Suppose  further  that  x  =  (0, 0, 0)  and  N(x)  =  (0, 0, 1).  Then  there  exists  a  neighborhood 
D  a(x)  E  S  of  x  such  that  Da(x)  admits  a  parametrization  given  by  the  formula 

Da(x)  =  {(«,  v,  g)  6  S|p(u,  v)  =  •  u2  +  k2  ■  v 2)  +  R(u,  u)},  (9) 

with  ki,  k-i  the  principal  curvatures  of  S  at  x,  and  the  function  R  satisfying  the  conditions 

|R(u,v)|  <  C  •  (u2  +u2)3/2,  (10) 

\Ru(u,v)\,  |J?„(u,u)|  <  C  ■  (u2  +v2),  (11) 

for  some  C  >  0  and  all  u2  +  v2  <  a2. 

3  Jump  Relations  of  the  Quadruple  Layer  Potential 

In  this  section,  we  will  derive  the  jump  relations  of  the  quadruple  layer  potential  on  a  regular 
surface  in  three  dimensions.  In  other  words,  for  x  €  S  let 

u(x  +  e  •  N(x))  =  Js  —  ^  (12) 

and  we  would  like  to  compute  the  limits  lime_>0+  u(x  +  e  •  N(x))  and  lim,-.  >0  u(x  +  e  •  N{x)). 
We  first  consider  the  simple  case  when  S  is  an  open  region  in  R2 . 
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3.1  Planar  Surface 


In  this  case,  we  may  assume  that  x  —  (0, 0, 0),  N(x)  =  (0, 0, 1),  t  =  (u,  v,  0).  Substituting  these 
equations  into  (3),  we  obtain 

d2G(x  +  e-N(x),t)  3e2 _ 1  (  . 

dN(t)2  (u2  +  V2  +  e2)5/2  (u2+v2  +  e2)3/2> 

and  thus 

.  ...  ..  f  d2G(x  +  c  N(x),t)  , 

u{x  +  e-  N(x))  =  J - - a(0  '  dt 

=  Is  {  (u2  +  V2  +  e2)5/2  ~  ( U 2  +  v2  +  e2)3/2  }  '  ^  ^  '  dU  '  dv' 

(14) 

The  following  lemma  states  that  the  quadruple  layer  potential  is  continuous  across  S  when  S 
is  flat. 

Lemma  3.1  Suppose  that  a  :  S  —t  M  is  twice  continuously  differentiable.  Then  the  function  u 
defined  by  (14)  satisfies  the  following  relation: 

lim  «(*  +  £  •  N(x))  =  Mm  £  {  ,„2  - 1^2)5/2  -  („2  +  J+  e2)3/2  }  ■  o(u,v)  ■  du  - iv 

=  -f'P'  I  („2  +(,2)3/2  ' "t”' «)  '  dU  ■  * 

=  -  ft  {  fs_B'  (^2)572  •  <2(«.  »)  •  *•  •  dv  ~  ^  ■  <2(0,  0)  }  , 

(15) 

with  De  =  {(u,t»)|u2  +  u2  <  e2}. 

Proof.  Obviously  the  problem  is  local,  and  we  may  assume  that  S  =  Dr  for  some  R  >  0. 

(a)  If  a(u,v)  =  1,  then  using  polar  coordinates  we  have 

U{x+e-  N(x))  =  JDr{  {u2  +  v2+  e2)5/2  -  (U2  +  y2  +  £2)3/2  }  '  1  '  ^  ^ 

2tt-R2 

-  (Ji;2  +  e2)3/2> 
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(16) 


and  thus 


Also 


limu(:r  +  e  •  N(x ))  = 
e— M)  V  V  "  R 


(17) 


-f.p.  I 

Js 


■  1  •  du  ■  dv 


s  ( U 2  +  V2)3/2 

=  -  lim  {  f  — — 1  ,  •  1  •  du  ■  dv  -  —  ) 

e~»0  \  Jdr-d t  (u2  +  v2)3/2  e  J 


(18) 


27T 

=  ~R' 

And  the  lemma  clearly  holds  in  this  case. 

(b)  Similarly,  the  lemma  holds  for  a  =  u  and  a  =  v  since  all  related  integrals  are  equal  to 
zero  by  symmetry  and  <r(0, 0)  =  0. 

(c)  We  now  consider  the  general  case  a  6  C2(S).  By  Taylor’s  theorem  we  have 

a(u,v)  =  cr(0, 0)  +  <ru(0,0)  •  u  +  cr„(0,0)  •  v  +  Ri(u,v),  (19) 

where  R]  satisfies  the  condition  that  for  all  u2  +  v2  <  82 

\Ri(u,v)\  <  M  ■  (u2  +  v2),  (20) 

with  J,  M  positive  real  numbers.  By  parts  (a)  and  (b),  we  only  need  to  prove  that  the 


lemma  holds  for  a  —  R\.  Since 

3e2 


( u 2  +  v2  +  e2)5/2  ( u 2  +  v2  +  e2)3/2 


|  •  Ri(u,v) 


<C- 


(u2  +  U2)l/2 

(21) 


for  all  e  <  eo  and  u2  +  v2  <  82,  the  integrands  in  (15)  are  only  weakly  singular  (and  hence 
absolutely  integrable)  even  if  e  =  0.  By  Lebesgue’s  dominated  convergence  theorem,  the 
order  of  limit  and  integration  can  be  interchanged  and  the  finite  part  integral  is  actually 
an  ordinary  (Lebesgue)  integral.  Hence,  the  lemma  is  proved. 

□ 

Next,  we  will  study  the  local  properties  of  the  quadruple  potential  near  a  regular  surface. 
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3.2  Local  Properties  of  the  Quadruple  Potential  near  a  Regular  surface 

Prom  now  on,  we  use  the  parametrization  given  in  Lemma  2.1,  i.e.,  x  —  (0, 0, 0),  N(x)  =  (0, 0, 1) 
and  for  t  €  Da(x )  C  S,  it  has  the  coordinates 


t  =  ( u,v,g(u,v )) 

1  (22) 
=  (u,  v, -(ki  ■  u2  +  k2  ■  v2)  +  R(u,  v)), 

with  k i,  k2  the  principal  curvatures  at  the  point  x  =  (0,0,0),  and  the  function  R  satisfying 
the  conditions  (10)  and  (11).  We  first  compute  the  quadruple  potential  in  the  neighborhood  of 
Da{x). 

Lemma  3.2  For  t  €  Da(x )  C  S, 

d2G(x  +  e-N(x),t)  _  3 (N(t),x  +  cN(x)-t)2 _ 1 

dN(t)2  \x  +  e  •  N(x)  —  f|5  \x  +  e  •  N(x)  —  t|3 

=  3  •  p2(u,v,e)  _  1  (23) 

J2(«,  v)  ■  d\(u,  v,  e)  di(u,v,e)’ 

where  the  functions  p\,  J,  d\  are  defined  by  the  formulae 

Pi(u,v,  e)  =  i(fci -u2  +  k2- v2)  +  Ru(u,v)  -u  +  Rv(u,v)  ■  v  -  R(u,v )  +e,  (24) 

J(u,  v)  =  \J\  +  k2  ■  u2  +  k2  -v2  +  2k\  -u  ■  Ru  +  2k2  ■  v  •  Rv  +  R2  +  R2,  (25) 

d\(u,  v,  e)  =  \x  +  e  •  N(x)  —  t| 

r.  TIT..  7“  _  ,  7  (26) 


+  V2  +  (l(fc,  ■  u2  +  ■  V2)  +  R( u,  v)  -  e)  . 


Proof.  The  first  equality  of  (23)  directly  follows  from  the  definition  of  the  quadruple  potential 
(3).  And  the  second  equality  follows  from  direct  computation  and  the  details  are  as  follows: 

x  +  e  •  N(x )  - 1  =  (-u,  —v, -\{k\  ■  u2  +  k2  ■  v 2)  -  R(u,  v)  +  e),  (27) 

tu  =  (l,0,k\  ■  u  + Ru),  tv  =  (0, 1,&2  -v  +  Rv),  (28) 

-  tu  x  tv  _  -  u-Ruiu,  v ),  —k2  ■  v  —  Ryiu,  v),  1) 


|  tu  X  tv  | 


J(u,  V ) 


(30) 


We  now  introduce  the  notation 


□ 


d(u,  v,e)  =  \/u2  +  n2  +  e2,  (31) 

The  following  lemma  provides  the  estimates  for  p\ ,  J,  and  di  to  be  used  in  the  proof  of  Lemma 
3.4. 


Lemma  3.3  there  exist  real  positive  numbers  C\,  C2,  C,  6,  to  such  that  for  all  t  G  D$(x)  and 
0  <  e  <  eo, 


|pi(u,u,e)2|  <  C  -d(u,v,t)2, 

\p\{u,  v,  e)2  —  e2  —  e  ■  ( ki  -u2  +  k 2  •  n2)|  <  C  ■  d(u,  v,  e)4, 
|  J(u,v)  —  1|  <  C  ■  d(u,v,e)2, 

C 1  •  d(u,v,e)  <  \di(u,v,e)\  <  C2  ■  d(u,v,e), 


1 

1 

3  •  e  •  (ki  ■  u2  +  k2  ■  n2) 

<c- 

1 

di  (u,v,  e)3 

d(u,v, 

e)3  2  -d(u,v,e)5 

d(u,  v,e)’ 

1 

1 

5  •  e  •  (fci  •  u2  +  k2  •  v2) 

<c- 

1 

d i(u,  v,  e)5 

d(u,v, 

e)5  2  -d(u,v,e)7 

d(u,v,e)3 

Proof.  The  estimates  (32)-(35)  directly  follow  from  the  definitions  of  p\,  J,  d\,  d 
condition  (10)  satisfied  by  the  function  R.  To  prove  (36),  we  first  note  that 

| d(u,  v,  h )  —  d\(u,  v,  fi)|  <C  ■  d(u,  v,  e)2, 

|d(-u,  v,  h)2  —  d\(u,  v,  h)2  —  t  •  (k\  •  u2  +  ^2  •  v2)\ 

=  2e  •  R(u,v)  —  f^(^i  •  u2  +  k2  •  v2)  +  R{u,  n)^ 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 
and  the 


(38) 

(39) 


<C  ■  d(u,v,  e)4. 
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For  the  sake  of  readability,  we  will  sometimes  denote  d\  (u,  v,  e),  d(u,  v,  e)  simply  by  di,  d 
respectively.  We  have 


1  e  •  (ki  ■  u2  +  k-2  ■  v 2)  •  (d?  +  di-  d  +  d2) 
d3  ■  d\  •  (di  +  d) 

(d2  —  d?  +  e  •  (fci  •  u2  +  A?2  •  ^2))  •  (d2  +  di  •  d  +  d2) 


df  d3 


d3  •  d3  •  (di  +  d) 


<C- 


d(u,v,e )’ 

where  the  inequahty  follows  from  (35)  and  (39).  Similarly 

(d2  +  di  •  d  +  d2)  3 

d3  •  (di  +  d)  _  2  ■  d2 
(d  -  di)  •  (3d?  +  6d  •  d?  +  4d2  •  di  +  2d3) 

2d2  •  d3  •  (dj  +  d) 

d(u,v,e ) 

where  the  inequality  follows  from  (35)  and  (38).  Furthermore, 


e  •  (&i  •  u2  +  &2  •  v 2) 
d(u,v,  e)3 


<  C. 


(40) 


(41) 


(42) 


Combining  (40)-(42),  we  obtain  (36).  Finally,  we  note  that  (37)  can  be  proved  in  an  almost 
identical  manner.  □ 

The  following  lemma  singles  out  most  singular  terms  in  the  quadruple  potential. 


Lemma  3.4  there  exist  real  positive  numbers  C,  8,  eo  such  that  for  all  t  €  Dg(x)  and  0  <  e  < 

eo, 


d2G(x  +  cN(x),t)  3e2 


+ 


<C- 


dN(t )2  d(u,  v,  e)5  '  d(u,  v,ej3 

3  ■  e  ■  (hi  •  u2  +  ~  v2)  15  •  e3  •  (ki  ■  u2  +  ki  ■  v2) 

2  •  d(u,v,t)5 

1 


2  •  d(u,  v,  e)7 


d(u,v,e) 

Proof.  Combining  (32)-(34),  we  have 

pi(u,v,e)2  e2 +t- (k\-u2 +  k2-v2) 


J{u,v )2  •  di(u,v,e)5 


di(u,v,  e)5 


<C- 


(43) 


d(u,v,  e) ' 


(44) 
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Combining  (23)  and  (44),  we  obtain 


d2G{x  +  e  •  N(x),  t ) 
dN(t)2 


3(e2  +  e  •  (ki  ■  u2  +  •  v2))  1 

di(u,v,  e)5  di(u,v,e)3 


<C- 


1 

d(u,v,e) 

(45) 


Finally,  we  obtain  (43)  by  combining  (36),  (37),  and  (45).  □ 

We  now  introduce  the  notation 

r(x ,  t)  =  r(u,  v )  =  \/  u2  +v 2.  (46) 


The  following  corollary  is  obtained  by  letting  e  =  0  in  (43). 


Corollary  3.5  there  exist  real  positive  numbers  C,  8  such  that  for  all  t  E  Dg(x), 

*)  .  1 

dN(t)2  r(x,t)3  ~  r(x,t)' 


(47) 


3.3  General  Regular  Surface 

We  now  consider  the  case  of  a  general  regular  surface.  Obviously  since  l/r(x,t)  in  (47)  is 
weakly  singular  (and  thus  absolutely  integrable),  we  obtain  the  following  lemma  by  Lebesgue’s 
dominated  convergence  theorem  and  Corollary  3.5. 


Lemma  3.6  Suppose  that  a  :  S 

d2G(x,  t) 
IS 


is  twice  continuously  differentiable.  Then 


f'P'^  dN(t)2 


a(t)  ■  dt  =  — f.p 


+ 


/,( 


’  Is  r3(*» 

d2G(x,  t ) 


*) 


a(t)  •  dt 


+ 


(48) 


a(t)  ■  dt , 


dN(t )2  r3(x,t)/ 

where  the  last  integral  is  interpreted  in  the  sense  of  an  ordinary  (Lebesgue)  integral. 


The  following  lemma  states  that  the  last  several  terms  in  (43)  are  all  “approximations  to  the 
identity”  (or  more  precisely,  a  constant) . 


Lemma  3.7  Suppose  that  f  €  Lp{Da{x)),  1  <  p  <  oo.  Suppose  further  that  x  belongs  to  the 
Lebesgue  set  of  f  (see,  for  example,  [5]  Chapter  I).  Then 


lim  f 

e^°+  Jd, 


e  •  u 


■  f(t)  -dt=^--  f(x), 


^o+JDa(x)d(u,v,e)5  •'v"/  3 


(49) 
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L 

im  [ 

■+°+  Jd, 


_  .  J'v  ■f(t)-dt=^--f(x), 

e~*°  +  Jdu(x)  d(u,v,e)5  3 


e3  •  u2 


2ir 


7  ‘  /(*)  •  dt  =  jg  '  /(*)> 


(50) 

(51) 

(52) 


e^°+  JDa(x)  d(u,  v,  e) 

r  p3  .  ..2  o„ 

lim  /  *  •  /(t)  -dt  =  —  -  fix). 

e^0+JDa(x)d(u,v,t)  !5 

(x  belongs  to  the  Lebesgue  set  of  f  if  lim^o  pr  Jdc(x)  1/(0  —  /(®)|  •  dt  =  0.  In  particular,  x  is 
a  point  of  the  Lebesgue  set  of  f  if  f  is  continuous  at  x.) 

Proof.  The  proofs  of  (49)-(52)  are  almost  identical  to  each  other  and  thus  we  will  only  prove 
(52).  Representing  dt  in  local  coordinates,  we  have 


f  d(  V  w  •/(*)•<**=  / 

JDa(x)  d{u,  v ,  e)7  y„' 


e3  •  v2 


2+t>2<a2  («2  +  +  h2)7/2 


■  f(u,  v )  ■  J(u,  v )  ■  dudv, 
(53) 


where  the  function  J  is  given  in  (25).  Obviously,  J  is  a  bounded  continuous  function  in 
u2  +  v2  <  a2  and  J(0, 0)  =  1.  Hence,  /( 0,0)  •  J(0, 0)  =  /(0,0)  =  f(x).  By  Theorem  1.25  of 
Chapter  I  in  [5],  we  only  need  to  prove  that 

v2  j  2tt 

—r - jr - -=nz  •  du-  dv  =  — . 

/K2  (u2  +  v2  +  l)7/2  15 


/ 

./R2 


(54) 


But  the  above  identity  can  easily  be  verified  by  evaluating  the  integral  via  polar  coordinates. 
Hence,  the  lemma  follows.  □ 

We  Eire  now  in  a  position  to  present  the  principal  result  of  the  paper. 

Theorem  3.8  Suppose  that  S  is  a  sufficiently  smooth  regular  surface  in  M3  and  that  cr  :  S'  — >  M 


is  twice  continuously  differentiable.  Then  the  quadruple  layer  potential  u 
the  formula 

d2G(x,  t ) 


u(x)  = 

Js 


a(t)  ■  dt 


is  dN(ty 

satisfies  the  following  jump  relations  at  x  £  S: 


defined  by 


(55) 
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(a) 


lim  u(x  +  e  •  N(x ))  =  lim 

e— ^0+  e-*-0+ 


L 


d2G{ x  +  e  •  N( x),  t ) 
dN(t)2 


•  a(t)  ■  dt 


2ir  ■  (h(x)  +  k2{x))  ■  u(x)  +  f.p.  jf  ■  <r(t) 


dt ; 


(56) 


(b) 


T  I  ,  AT/  w  v  f  d2G(x  -f  e  •  N(x),t)  .  . 

e-+o-  v  v  "  t-+o-Js  dN(t)2  w 

=  -2t r  •  (fci(as)  +  fc2(x))  •  a(x)  +  f.p.  £  '  CT(*)  '  dt> 

(57) 

wi/i  /cj  (a:),  &2(x)  the  principal  curvatures  of  S  at  x.  (Note:  k\(x)  +  k2(x)  is  the  trace  of 
the  second  fundamental  form  of  S  at  x,  which  changes  sign  ifN(x)  changes  sign  (i.e.,  S 
changes  its  orientation).) 

Proof.  The  proofs  of  (56)  and  (57)  are  almost  identical,  so  we  will  only  prove  (56).  By  lemma 
3.4,  we  have 


d2G(x  +  e  •  N(x),t ) 
dN{t)2 


•  cr(f)  •  dt  =  f  I  — rr  -  ~t> — - — vs}  •  <x(f) 

Js\d(u,v,e )5  d(u,v,e)3j 


L 

f  f  3  •  e  •  (fci  •  u2  +  k2  ■  v 2) 

+  Js  l  2  -d(u,v,e)5  +  2  •  d(u,v,  e)7 

+  /  R(u,v,e)  ■  cr(t)  ■  dt 

Js 


•  dt 


15  •  e3  •  ( k\  ■  u2  +  k2 


V)| 


cr(t)  •  dt 


:=  h  +I2  +  I3, 


(58) 


where  the  function  R  satisfies  the  conditions 


|i?(n,  v,  e)|  <  C  ■ 


1 

d(u,  v,  e) ’ 


R(u,  v,  0) 


d2G(x,t)  _ 1_ 

dN(t)2  +  r3(x,  t)' 


(59) 

(60) 
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We  now  analyze  the  above  three  items  separately.  First,  by  Lemma  3.1,  we  have 


-,P./ 

Js 

Second,  by  Lemma  3.6,  we  have 


lim  I\ 

e— lO-f 


r(x,ty 


a(t)  •  dt. 


(61) 


lim  Io 
£->0+ 


2n  ■  ( ki(x )  +  k2(x))  ■  a(x). 


(62) 


Third,  by  Lebesgue’s  dominated  convengence  theorem,  we  have 

.5?+  h  =  ls  R(u'  0  )■'«■)■*  =  fs{  } '  "M '  *•  (63> 

Finally,  we  obtain  (56)  by  combining  (48),  (58),  (61)-(63).  □ 
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